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Abstract 

Analytical methods are used to prove the existence of a periodic, symmetric 
solution with singularities in the planar 4-body problem. A numerical calculation 

Qand simulation are used to generate the orbit. The analytical method easily extends 
to any even number of bodies. Multiple simultaneous binary collisions are a key 
^ <~i feature of the orbits generated. 
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^ 1 Introduction 

The n-body problem of celestial mechanics is one of the most important problems in 
the field of dynanaical systems. The following differential equation 

O pi \?i-?iV 



> 



gives a mathematical description of the planar n-body problem, where p, £ K denotes 
the position of the rth body having mass m,-. All derivatives are taken with respect to 
time t. The potential energy of the system is given by 

1<(<7<«IP« Pjl 

and the kinetic energy is given by 

T =\j;^mi\pi\\ (3) 

^ (=1 
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Linearly stable symmetric periodic orbits are one aspect of the n-body problem. 
The elliptic Lagrangian triangular periodic orbits are linearly stable for certain values 
of eccentricity and the three masses lISl. lTTl . The Montgomery -Chenciner figure-eight 
orbit for three equal masses ||3l, ^ has been shown by Roberts [ lOJ to be linearly sta- 
ble by an innovative symmetry reduction technique he developed. 

Singularities are another particular aspect of the n-body problem. Binary collisions, 
triple collisions, etc,, are discussed at length in |13|. The Simultaneous Binary Colli- 
sion (SBC) problem has been widely studied as well, both analytically and numerically. 
Simo Q showed that the block regularization in the cases of the n-body problem which 
reduce to one-dimensional problems is differentiable, but the map passing from initial 
to final conditions (in suitable choices of transversal sections) is exactly C**/^. Ouyang 
and Yan lfT2ll give another approach for the regularization and analyze some properties 
of SBC solutions in the collinear four-body problem. Elbialy Q studied the nature of 
the collision-ejection orbits associated to SBC. 

Schubart |14] combined these two aspects to produce a singular linearly stable pe- 
riodic orbit in the three-body equal mass collinear problem. The motion of the middle 
mass regularly alternates between binary collisions with each of the outer two masses. 
His work was subsequently extended to the unequal mass case by both Henon |5 1 and 
Hietarinta and Mikkola |6|. Sweatman |15| later extended this work to a four-body 
periodic solution in one dimension, with bodies alternating between SBC of the outer 
mass pairs and binary collision of the inner two masses. 

In this paper, we present the analytic existence of a family of singular symmet- 
ric periodic planar orbits in the four-body equal mass problem. The initial conditions 
of these orbits are symmetric in both positions and velocities, which lead to periodic 
simultaneous binary collisions with each of the four masses alternating between col- 
lisions with its two nearest neighbors. Due to the abundance of symmetries present 
in the initial conditions, we can reduce the number of variables needed to just four - 
two for representing position and two for representing momentum. In contrast to its 
one-dimensional counterparts, the proof for existence of this orbit is surprisingly sim- 
ple. We begin in Sect ion | 2.1| by giving a description of the proposed orbit and prove its 
existence. In Section |2.2| we present the numerical methods used to produce the initial 
conditions that will lead to this orbit. In Section |3] we consider variants on this orbit, 
giving a family of orbits with singularities for an even number of equal masses. 

Since the initial submission of this paper, we have been doing additional work with 
Dr. Lennard Bakker (Brigham Young University) and Dr. Gareth Roberts (College of 
the Holy Cross) implementing Robert's linear stability technique as presented in lITOl . 
After precisely defining the symmetries that are present in the regularized coordinates, 
it is shown that the group of symmetries in the orbit is isomorphic to the dihedral group 
D4. Further, as a consequence of Robert's technique, we have shown that the four-body 
planar orbit presented in this paper is linearly stable ||T|. Further work has also been 
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done on orbits in this family with ahernating unequal masses (T\. Rather than a single 
mass parameter, the bodies have masses mi, 1112, mi, m2 as numbered moving counter- 
clockwise through the plane. Since some symmetries have been lost by this change in 
masses, it is necessary to choose two initial condition parameters as well as two initial 
velocities. Although numerically this is not a difficult problem, an analytical technique 
will require much more work. 



2 The Proposed Orbit 
2.1 Analytical Description 



T 
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Figure 1: On the left, we illustrate the initial conditions leading to the 
four-body two-dimensional periodic SBC oribt. On the right, the orbit is 
shown. 



Initially we focused on finding a symmetric, periodic SBC orbit for four equal 
masses in two dimensions. Without loss of generality, we assume that the orbit be- 
gins with the four bodies lying at (± 1 , 0) and (0, ± 1 ) with initial velocities (0, ±v) and 
(±v, 0), respectively, where v £ (0,+°°). For convenience throughout the rest of the 
paper, we number the bodies 1 to 4 as in Figure 1 . 

The singularity of SBC in this problem is not essential. For a better understand- 
ing of the behavior of the motion of the bodies in a neighborhood of a collision, the 
standard technique is to make a change of coordinates and rescale time. In the new 
coordinates, the orbits which approach collision can be extended across the collision 
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in a smooth manner with respect to the new time variable. This technique is called 
regularization. In our problem, the regularization describes the behavior of the bodies 
approaching and escaping collisions, similar to the collisions of billiard balls. 

Due to the symmetry of the initial conditions and the equations governing the mo- 
tion of the bodies, the symmetry that is present in the initial conditions is maintained 
in the regularized sense. 

Main Theorem. Let E = T be the total energy and m be the mass for each of the 
four bodies. For any E <0 and m > 0, there exists a symmetric, periodic, four-body 
orbit with SBC in M?. 

Without loss of generality, we can assume m = 1 and the initial positions are as 
illustrated in Figure 1. The proof will be given at the end of this section. 

Let to be the time of first SBC. For f £ [0,fo), let the coordinate of body 1 be 
(xi,X2). By symmetry, the coordinates of bodies 2, 3, and 4 are {x2,xi), (— xi,— X2) 
and {—X2,—xi), respectively. Using equation ([Tji, the acceleration of a body at point 
{xi,X2) is given by: 



{xi,X2) 



(xi -jt:2,Jt:2 -xi) ^ (2x1,2x2) ^ (xi +X2,xi +X2) 
(2(xi-X2)2)i (4x2+4x2)1 (2(xi+X2)2)i 



(4) 



^1 -^-^-^2) 

We now perform the regularization of the system. The system has the Hamiltonian: 

1/ 2 2> V2 V2 1 

H=-iw^,+wl)-^ ^ ^= (5) 

8 XI -X2 X1+X2 ^jcl+xj 

where wi = 4xi and W2 = 4x2 are the conjugate momenta to xi and X2. Note that SBC 
happens when xi — ±X2. We introduce a new set of coordinates: 

^i=xi-X2, q2^xi+X2. 

Their conjugate momenta pi are given by using a generating function F — {xi~X2)pi + 

(xi +X2)p2- 

W\^P\+P2, W2 = P2-Pl- 

The Hamiltonian corresponding to the new coordinate system is 

1, 7 7s \/2 \/2 \/2 
H=-Ap\ + pI)- — -- (6) 

Following the work of Sweatman ifTSl . we introduce another canonical transforma- 
tion: 

qi^Q}, Pi^2QiPi (/=1,2) 
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with Qi > 0. We also introduce a new time variable s, which satisfies ^ = qiqi- This 
produces a regularized Hamiltonian in extended phase space: 

as 

where E is the total energy of the Hamiltonian H. 

The regularized Hamiltorrian gives the following differential equations of motion: 



G'l - 7,PiQl (8) 



Q'l = IpiQI (9) 



1 7 r- 2V2QiQ^ 1^Q\Q 

p; = — pIQ\ +2V2Q1 + ; ^ ^- ^' ■ 



-i---^2^i^^v^^i^^^=- ^ '-'"f^ +2EQiQl (10) 



with initial conditions 

ei(0) = l, 62(0)-!, Pi(0) = -4y, P2(0)=4y (12) 
where derivatives are with respect to s, and E is the total energy of the Hamiltonian H. 



Theorem 1. Let sq be the time of the first SBC in the regularized system. Then sq is a 
continuous function with respect to the initial velocity v. Furthermore, 

P2[t0) - 



262(^0,0 

is also continuous with respect to v. 



Proof. At the first SBC, Qi (sq) = 0, and 22(50) = ^/qi = V-^i +X2 > 0. Our goal is 
to show that p2{to) is a continuous function with respect to v. 

Because r = at i = io, Pi {sq) = -4</2 from Since F is regularized, the 
solution Pi — Pi{s, v) and Qi — Qi{s,v) are continuous functions with respect to the two 
variables s and v. At time s — sq, 

o = ei(5o(v),v). 
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To apply the implicit function theorem, we need to show that 

wn 

-^(io,v) ^0. 

Fmm^ 

^(^o,v) = ^-PiQl I (,„,,„)= -\</2Q2{sof < 0. 

So SQ = so(v) is a continuous function of v. Therefore both P2{sq,v) and Q2{so,v) are 
continuous functions of v. Further, since Q2{so,v) > 0, p2(to) is also a continuous 
function of v. □ 

Theorem 2. There exists a v — vq such that x\ (fo) +i2(fo) = j/'2(fo) 0, where to is 
the time of the first SBC, i.e. the net momentum of bodies 1 and 2 at the first SBC is 0. 

The outline of this proof is as follows: We will show that there exist vi and V2 such 
that x\ +X2 is negative at SBC for v = vj and positive at SBC for v = V2. The result 
then follows by Theorem [T| 



Proof. Consider Newton's equation before the time of the first SBC: 

.. _ X2-XI 2x1 XI +X2 

'""2V2(x,-X2)3 8(x2+x2)3/2 2V2(xi+X2)3' ^ ^ 

.. _ Xi -X2 2X2 Xi +X2 

'''"2V2(xi-X2)3 8(x2+x2)V2 2V2(xi+X2)3' 

Therefore, 

Xl +X2 1 

X'l +X7 — ^TTT^ 1= T < 0, (15) 

- 4(x2+x2)3/2 V2(X1+X2)2 

which means xi +X2 is decreasing with respect to t. 

At the initial time f = 0, xi = 1, X2 = 0, xi = 0, and X2 = v. Note that for v e (0, 0°), 
there is no triple collision or total collision for f G [0,fo], where to is the time of the first 
SBC, as a triple collision implies total collapse by symmetry. Also, from the initial 
time to fo, < X2 < Xl < 1, < xi +X2 < 2, and Xj +X2 < 4. 

Let ^(f) — Xl (f) +X2(f). Then for any choice of v, y{t) < and < ^(f) < 2 hold 
for any t E [Q,to]. In other words, y{t) is decreasing with respect to t. 

First, we will show that there exists vi such that y{to) < 0. When v = the four 
bodies form a central configuration and, as a consequence, the motion of the four bod- 
ies leads to total collapse. Consider the time interval t E [0,fo/2). In this interval, the 



differential equations ( 13 1 and ( 14i have no singularity, and y{to/2) < 0. By contin- 
uous dependence on initial conditions, y{to/2) —xi{to/2) +X2(fo/2) is a continuous 
function with respect to the initial velocity v. When v = 0, xi (to/2) < 0, X2(ro/2) = 0, 
which gives y{tQ/2) < 0. Therefore, there exists a 5 > 0, such that y{to/2) < holds 
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for any v € (—5,5). 



Choose vi — 5/2, then j(fo/2) < 0. Because y{t) is decreasing with respect to t, 
j(fo)<j(fo/2)<0. 

Next we will show that there exists V2 such that ^(ro) > 0. Note that as v 0°, 
lim3;(fo) = lim xi (fo) +X2(fo) =2 

and 

lim y{t()) = 00. 

Therefore there exists some positive value V2, such that y{to) > 0. □ 

Proof of the Main Theorem. From Theorem[2] we know there exists an initial velocity 
V = vo such that ii (fo) +i2(fo) = 0. Let {Pi,P2tQ\ , 62} for s e [0,io] be the solution 
in the regularized system corresponding to the orbit from t = to t = to. Following 
collision, consider the behavior of the first and second bodies. Assume their velocity 
was reflected about the y = x line in the plane. In the new coordinate system, this 
corresponds to a new set of functions 

{ -Pi (Iso -s),-P2{2sQ-s),-Qi {2so -s),-Q2 (2io -s)} 

for s G [so,2so]. We can easily check that 

{-Pi (Iso - s),~P2{2sQ - s),-Qi {2so - .9), -62(2^0 - s)} 



for s e [io,2io] is also a set of solutions for equations ([8]) through (111 with initial con 



ditions at s = sq. Also, {Pi(i),P2(^),2i(.s'),22('^)} for s G [io,2io] satisfies equations 



([8J1 through (111 with the same initial conditions at s — sq. Note that equations (|8]l 
through ( [TT| with initial conditions at s = sq have a unique solution for any choice of 
V e (0,00). Then by uniqueness, the orbit for s E [sq,2so] must be the same as the orbit 
for s e [Q,S(,] in reverse, i.e. 

P,{s) = -P,(2.9o - s),Qi{s) = -Q,{2sa - s) 

for s e [0,io]- Therefore at time s — 2sq, bodies 1 and 2 will have returned to their ini- 
tial positions with velocities (0, — v) and (— v,0) respectively. Similarly, at time s = 2sq, 
bodies 3 and 4 will have also returned to their initial positions with velocities (0, v) and 
(v,0) respectively. 



Next, we use symmetry and uniqueness to show the orbit from s = 2sq to s = 4so 
and the orbit from s — Q Xo s — 2.so will be symmetric with respect to the y— axis. 
Compare the motion of body 2 and body 3 from s = 2s() to i = 4^0 with the motion 
of body 2 and body 1 from time s = to s = 2so. The initial conditions of body 3 at 
s ~ 2so and the initial conditions of body 1 at i = are symmetric with respect to the 
y— axis. Also the initial conditions of body 2 at s = 2so and the initial conditions of 
body 4 at .? = are symmetric with respect to the x— axis. Therefore, by uniqueness. 
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the orbit of bodies 2 and 3 from s = 2so to i = 4so and the orbit of bodies 1 and 2 from 
s = to s = 2so must be symmetric with respect to axis. Therefore, the orbit of 
bodies 1 and 4 from s ~ 2so to i = 4so and and the orbit of bodies 3 and 4 from s ^0 
to s = 2sq are symmetric with respect to the y— axis. Hence, at 5 = 4^0, the positions 
and velocities of the four bodies are exactly the same as at 5 = 0. Therefore, the orbit 
is periodic with period s = 4so. □ 

It is worth noting here that the previous proof implies a time-reversing symme- 
try for the periodic orbit. This provides further evidence for the conjecture made by 
Roberts fTOl, stating that linearly stable periodic orbits in the equal mass n-body prob- 
lem must have a time-reversing symmetry. (Linear stability of this orbit is shown in 

m.) 

2.2 Numerical Method 

As we are searching for a periodic orbit of the «-body problem, we assume the value 
of the Hamiltonian needs to be negative. Using the initial positions of the four bodies 
described earlier, it is not hard to find the potential energy at f = 0: 

U = 2V2 + l. 

Then, acting under the negative Hamiltonian assumption: 



Since all masses are equal, if we require that the velocities of each body are equal in 
magnitude, we obtain: 



with v^ax defined to be the value of v such that the value of the Hamiltonian is zero. 
Define 9 — This parameter is used in the numerical algorithm. 



At this point it becomes necessary to find out just how much kinetic energy is re- 
quired to obtain the periodic orbit. Since we know suitable bounds on the velocity 
parameter (0 S (0, 1)), we can search the interval numerically. We use an n-body sim- 
ulator with the initial positions previously described. The simulation is ran until one 
SBC occurs. For simplicity, we consider only the collision between the first and second 
bodies in the first quadrant. Summing their velocities immediately before the collision 
gives a vector running along the line y — x (due to symmetry), with both components 
having the same sign. The magnitude of this vector is given in Figure 2. Negative 
magnitudes represent vectors with both components less than zero. 



2%/2+l > ^ 



i=l 



n 





(16) 
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Figure 2: The magnitude of the net velocity of the first two bodies (ver- 
tical axis) at the time of collision for various values of 9 (horizontal axis) 



Next, a standard bisection method is used to find the amount of energy required 
to cause the net velocity at collision to be zero. Using the initial interval 9 G [0, 1] 
and iterating to a tolerance of 10^'^, the correct value of 9 was found to be 9 = 
0.46449539554694. 

It is worth noting that both the proof of existence and the numerical method do 
not guarantee the uniqueness of this orbit. Numerical simulations demonstrate that for 
values of 9 near the correct value, the orbit remains for a significant length of time with 
the paths of the bodies lying in a "fattened" annular region roughly the shape of the 
original orbit. Near the extreme ends, the orbit experiences near total-collapse and fall 
apart rapidly. Although we do not focus on these questions just yet, a more thorough 
study of the dynamics could prove to be quite interesting. 

3 Variants 

3.1 Orbits of more than four bodies 

The same technique can be adopted to find similar orbits for any arbitrary even number 
n. A key feature of these orbits will be higher numbers of simultaneous binary colli- 
sions. For a given value of «, initial positions are given by spacing the bodies evenly 
about the unit circle. The potential energy (and the value of Vmax) is found numeri- 
cally by iterating over each pair of planets and summing the reciprocal of the distances 
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between them. (Recall that all m, = 1.) Velocities are then assigned to the bodies in 
alternating counter-clockwise and clockwise directions, initially tangent to the circle. 
Again we consider the collision between the first and the second bodies. Although the 
net velocity of the two at collision will not lie along the y — x line, the components of 
this vector will both have the same sign. The magnitudes of the net velocity between 
the first two bodies at initial collision are shown in Figure 3 for various values of n. 
Lower curves in the graph correspond to higher values of n. Again, negative magni- 
tudes correspond to both components being negative. 




Figure 3: Curves showing the magnitude of the net velocity of the first 
two bodies (vertical axis) at the time of colhsion for various values of 
(horizontal axis) for n ~ 4,6,8, 10, 12. 



Pictures of the orbit for n = 6 and n = 8 are shown in Figure 4. It is readily seen 
that as n increases, the shape of the orbit more closely approximates a circle. 
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Figure 4: The six- and eight-body two-dimensional periodic SBC orbits. 
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